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Parameters of nuclear density distributions are derived from least-squares fits to strong interaction 
observables in exotic atoms. Global analyses of antiprotonic and pionic atoms show reasonably good 
agreement between the two types of probes regarding the average behaviour of root-mean-square 
radii of the neutron distributions. Apparent conflict regarding the shape of the neutron distribution 
is attributed to different radial sensitivities of these two probes. 
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I. INTRODUCTION 

The density distribution of protons in nuclei is considered known as it is obtained from the nuclear charge distribution 

by unfolding the finite size of the charge of the proton. The neutron distributions are, however, generally not known 

to sufficient accuracy. A host of different methods has been applied in studies of root-mean-square (rms) radii of 

neutron distributions in nuclei but the results are sometimes conflicting, see e.g. In the present work we focus 

^ ' on antiprotonic and on pionic atoms as a source of information on neutron densities. We deduce average properties 
OO ' 

^-j- I with the help of global analyses of about 100 data points in each case, covering the whole of the periodic table. 
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Reasonably good agreement is obtained between the two types of exotic atoms when considering rms radii of the 
neutron distributions. A conflict regarding the shape of the distributions is most likely due to the different radial 
sensitivities of the two probes. 



II. METHOD 

Strong interaction level shifts and widths in exotic atoms, formed by the capture of a negatively charged hadron 
into an atomic orbit, are calculated with the help of an optical potential inserted into the appropriate wave equation. 
The simplest class of optical potentials V op t is the generic tp(r) potential, which for underlying s-wave hadron-nucleon 
interactions assumes the form: 

2pV opt (r) = -4tt(1+ ^-J?){Mpn{r) + Pp (r)] +T z b 1 [p n (r) - Pp (r)]} . (1) 

Here, p n and p p are the neutron and proton density distributions normalized to the number of neutrons N and number 
of protons Z , respectively, M is the mass of the nucleon and r z = +1 for the negatively charged hadrons considered 
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in the present work. When handling many different nuclei over the periodic table it is necessary to represent the 
densities by approximate distributions, usually chosen as the the two-parameter Fermi distribution (2pF). With the 
proton densities considered known, we focuse on the differences between the neutron and the proton distributions. 
A linear dependence of r n — r p , the difference between the rms radii, on (N — Z)/A has been employed in p studies 
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namely 



N - Z 

r n ~r p = 7 — -r h S , (2) 

with 7 close to 1.0 fm and 5 close to zero. This parameterization is adopted here. In order to allow for possible 
differences in the shape of the neutron distribution, the 'skin' and 'halo' forms of Ref. Q were used, as well as an 
average between the two. We adopt a 2pF distribution both for the proton (unfolded from the charge distribution) 
and for the neutron density distributions 

I \ P0n,0p / \ 

PnA r ) = T~i Ti u — w V ■ ( 3 ) 

1 + cxp((r - Rn,p)/a n ,p) 

Then for each value of r„ — r p in the 'skin' form the same diffuscness parameter for protons and neutrons, a n = a p , 
is used and the R n parameter is determined from the rms radius r n . In the 'halo' form the same radius parameter, 
R n = R P , is assumed and a n is determined from r n . In the 'average' option the diffuseness parameter is set to be the 
average of the above two diffuseness parameters, a n vc = (a p + a^)/2, and the radius parameter R n is then determined 
from the rms radius r n . In this way we can test three shapes of the neutron distribution for each value of its rms 
radius all along the periodic table. These shapes provide sufficient difference in order to be tested in global fits. The 
results below are presented as the best fit \ 2 values vs. the neutron rms radius parameter 7. 

III. RESULTS 
A. Antiprotonic atoms 

Figure Q] shows results of global fits to 90 data points from measurements of X-rays by the PS209 collaboration 0|. 
It is seen that the 'skin' shape for the neutron density distribution is unfavoured and that a finite range for the pN 
interaction leads to significant improvements in the fits. Comparing the zero-range (ZR) results of the left-hand side 
of Fig[T] with the corresponding ZR results of Ref. [2| we note that the values of the best-fit parameter 7 are quite 
different. This is due to the use in Ref.Q] of fixed values for the complex parameter bo for all values of 7 whereas we 
re-fit these phenomenological parameters when 7 changes. In fact, with the fixed values for 60 taken from Ref. [(J] one 
obtains a value for 7 that merely represents an average over the neutron densities used in Ref. [6| to derive those fixed 
values. 

The best fit finite-range potential produces a \ 2 P er degree of freedom of about 2, which is most acceptable 
considering the simplicity of the model and the extent of the data. We also note that on the basis of values of \ 2 it is 
impossible to distinguish between the 'halo' and the 'average' shapes, which lead to somewhat different values of the 
rms radius parameter 7. Similar analyses of the radio-chemical data and of the combined X-rays and radiochemical 
data lead to very similar conclusions. It is interesting to note that the isovector parameter b\ turns out to be consistent 
with zero [4] . This is in full agreement with the result of analysing separately the radio-chemical data where we find 
that the best fit is obtained for a ratio of 0.99±0.07 for the absorption on a neutron to the absorption on a proton. 
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FIG. 1: Zero range (left) and finite range (right) fits to strong interaction shifts and widths in antiprotonic atoms. For the 
latter the best fit is obtained for a rms radius of the pN interaction of 1.1±0.1 fm. 



B. Pionic atoms 



Values of rms radii of neutron distributions for several nuclides had been derived from pionic atoms more than 
two decades ago Q|. The first to derive rms radii from extensive data sets of strong interaction observables in pionic 
atoms were Garcia- Recio et al. [7]. Using the 'skin' shape for the neutron densities they presented a list of rms radii 
which, when interpreted with the present formulation, leads to 7 = 1.06 ± 0.34 fm for their semi-theoretical model 
and 7 = 0.97 ± 0.12 fm for the phcnomcnological model of Meirav et al. [||. 

Figure [2] shows results of global fits to 100 data points for pionic atoms, using the latest pion-nucleus potential 
with energy and density dependence in the s-wave term and finite range in the p-wave term, see [{J for details. 
The left hand side shows that the 'skin' shape for the neutron density distributions yields the lowest \ 2 value. The 
horizontal band represents the value of x 2 per degree of freedom, which indicates the statistical significance of the 
fits and determines the uncertainty of the derived parameters, 7 in this case. On the right hand side we compare 
data for Pb with predictions made with the global parameters of the left hand side. The results are consistent with 
the global analyses but the uncertainties are obviously considerably larger. This is typical also of antiprotonic atoms, 
demonstrating the limited accuracy of a single element analysis. 



IV. DISCUSSION 



Table U summarizes the above results regarding the shapes of the neutron density distributions, within the simple 
2pF parameterization, and the values of the parameter 7 of Eq.@, which determines the dependence of the rms radius 
on the neutron excess parameter (N — Z)/A. The other parameter was held fixed at S = —0.035 fm. For antiprotonic 
atoms the two shapes of the neutron distributions lead to almost the same quality of fit and it is impossible to 
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FIG. 2: Left: global fits to pionic atom data. S, A and H represent the skin, average and halo shapes, respectively. The 
horizontal band represents \ 2 P er point, see text. Right: similar to the left but for pionic atoms of Pb. 



TABLE I: Summary of results showing number of points and best fit values of \ 2 
source shape of p n N \ 2 7 (f m ) 



'halo' 90 196 0.9±0.1 

'average' 90 198 1.25±0.15 



P 
P 

7T~ 



'skin' 



100 173 1.0±0.1 



distinguish between the two on the basis of the values of x 2 . The value of 7 from pionic atoms is consistent with 
either values obtained from antiprotonic atoms but there seems to be a conflict considering the shape of the neutron 
density. 

To look for the source of this conflict it is necessary to look into the radial sensitivities of the two probes, which 
have vastly different absorption cross sections in nuclear matter. This is done with the functional derivative method, 
introduced originally in connection with kaonic atoms 10] and later used also for pionic and antiprotonic atoms Q|. 
It is shown in [9| that pionic atom data are sensitive to nuclear densities around the 50% region of the central density 
whereas antiprotonic atom data depend on the density at the extreme periphery where the densities are well below 
10% of the central density. 

Figure [3] shows comparisons between the three versions of the 2pF neutron density and two, more physical, neutron 
densities for 208 Pb: (i)A single particle (SP) neutron density obtained by filling in of single particle levels in a common 
potential. (ii)A neutron density from a more sophisticated RMF calculation [llj]. All five densities have the same 
rms radius. On the left hand side is indicated the sensitive region for pionic atoms and on the right hand side is 
indicated the sensitive region for antiprotonic atoms. It is evident that in the pionic atoms region the skin and the 
average shapes are closest to the SP and the RMF densities whereas in the antiprotonic atoms region the skin shape 
deviates markedly from the two more physical models. At very large radii the halo shape is very close to the RMF 
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FIG. 3: SP and RMF neutron densities in 208 Pb compared with the three options of 2pF parameterizations (S for skin, H for 
halo and A for average). Left: covering the pionic atoms sensitive region. Right: covering the antiprotonic atoms sensitive 
region. All five densities have the same rms radius. 



density but at the antiprotonic atoms region of sensitivity the average shape is equally good. The conclusion is that 
the 2pF model is not necessarily able to reproduce more realistic densities over a broad range of density values. In 
contrast it seems that the values of rms radii are less sensitive to the model used. We conclude that there is no real 
conflict between antiprotonic and pionic atoms and that on the average the rms radii of neutron density distributions 
in nuclei may be represented by Eq.([2]) with 7=1.Q±0.1 fm and 5 = —0.035 fm. Applying this to 208 Pb which is a 



most studied nuclide, we find r n — r p =0.18±0.02 fm, in very good agreement with Ref. [12 1. 
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